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Abstract 

We investigate the long-time behavior of weak solutions to the thin-film type equation 

v t = (xv - vv xxx ) x , 



which arises in the Hele-Shaw problem. We estimate the rate of convergence of solutions to the 

24 I 



Smyth-Hill equilibrium solution, which has the form -kAC 2 — £ 2 )+, in the norm 



11/11™,!= /(l + N 2m )|/(x)| 2 dx+ / |/*(x)| a d*. 
JR JR 

We obtain exponential convergence in the ||| • | mj i norm for all m with 1 < m < 2, thus obtaining 
rates of convergence in norms measuring both smoothness and localization. The localization is 
the main novelty, and in fact, we show that there is a close connection between the localization 
bounds and the smoothness bounds: Convergence of second moments implies convergence in 
the H 1 Sobolev norm. We then use methods of optimal mass transportation to obtain the 
convergence of the required moments. We also use such methods to construct an appropriate 
class of weak solutions for which all of the estimates on which our convergence analysis depends 
may be rigorously derived. Though our main results on convergence can be stated without 
reference to optimal mass transportation, essential use of this theory is made throughout our 
analysis. 
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1 Introduction 

1.1 The primary Lyapunov functional 

The following one-dimensional fourth-order nonlinear degenerate parabolic equation 

m = - {uu xxx ) x , x£R,t>0, (1.1) 

with 

u(x,0) = u (x) > 0, iGR, (1.2) 

arises as the particular case of the thin-film equation in the the Hele-Shaw setting (21 [Qj, 110] . 

Equation (jl.ip has a well-known scale invariance and corresponding self-similar solutions [5|I12|: 
Given a solution u(x,t) of (jl.ip . define «(x, t) by 

u(x, t) := a(t)u(a(t)x, b(t)), (1.3) 

where a(i) = e* and 6(t) = (e 5t — l)/5. Then v(x, t) satisfies: 

v t = (xv - vv xxx ) x , x G R,i > 0, (1.4) 

v(x,0) = u {x), x£R. (1.5) 

Both equations (jl.ip and (jl.4p are conservation laws, so that the total mass is conserved. That 
is, in the case of (|1 .4|) . 

M = v(x,t)dx 
Jr 

is conserved. 

Also, both equations describe a gradient flow [5j: For fjl .41) . define the energy functional 

= o / ( -y ^( :r ) + ^ 2 ^(^)) dx. 
Then, (|1.4p can be written as 

-(•(TO." 

It follows from (jl.6p that for sufficiently regular solutions of (jl.ip . 



where 



5v 



dx . 



Thus, as long as the solution v is non- negative, E[v] will be monotone decreasing. 

The physical interpretation of u as the height of a thin film [2j[9l[T0] suggests that u, and hence 
v, must be non-negative for all physically meaningful solutions. However, there is no maximum 
principle to provide a-priori assurance of non- negativity for this fourth order equation, and it is an 
open question whether certain solutions, necessarily non-physical, can become negative. In what 
follows, we shall construct and analyze a particular class of non-negative solutions. 

In the rest of this introduction, we restrict our attention to non- negative solutions for which 
E[u] is monotone decreasing. 



CU May 25, 2010 



3 



It is easy to determine the minimizers of E[v], each of which is of course a steady state solution 
of (jl.4|) . and thus determines a self-similar solution of (jl.ip : These minimizing steady states [12] 
are the one parameter family of functions 

v^\x):=±(C 2 -x 2 ) 2 + , (1.7) 

where the subscripted + indicates the positive part, and the parameter C determines the total 
mass of i/ 00 -*. Because mass is conserved, we suppress C in our notation. 
Define the relative energy E\v\v^°°^\ by 

E[v\v^] = E[v] - E[v^} . (1.8) 

It turns out that this is a good Lyapunov functional for the steady states i>(°°): 

1.1 LEMMA. Let v have the same mass as v^°°\ Then: 

E[v\v^} > \ I \v x -v^\ 2 dx + If \x\ 2 vdx . (1.9) 

Proof: One finds, after one integration by parts, and using the fact that v and have the same 



mass. 



E\V\V {00) } -If IVr-vi^^dx 



2 



4 / \x\ 2 v^dx + \c 2 [ v^dx- [ (v^fdx 
2 Vm 6 ./» ./w 



{V°°)=0} 



-C 2 - -\x\ 2 ) 
6 2 1 



By direct computation, the three integrals on the second line cancel exactly, and on {u(°°) = 0}, 
\x\ 2 /2-C 2 /6 > \x\ 2 /3. □ 

Thus, i^l?/ 00 )] is monotone decreasing along the evolution described by (|1.4p . and by p. 91) . 
whenever ^[ulf (°°)] is small, v is close to v(°°\ which is what we mean by saying that v^ 00 - 1 ] is 
a good Lyapunov functional. Let us express this quantitatively. We shall use the following norms: 

1.2 DEFINITION. For any smooth compactly supported function / on R, and any m > 0, define 

I* ! = / (1 + \x\ 2m )\f(x)\ 2 dx + [ \f x (x)\ 2 dx . (1.10) 



We extend then the norm ||| • ||| mj i to its natural Hilbert space domain. 
1.3 LEMMA. Let v > have the same mass M as Then: 

Hoc < \m^{2E[v]) 1 ^ , (1.11) 
and there exist explicitly computable constants K± and Ki depending only on M such that 

\\v-v { - oo) l 2 ll <K l {E[v\v { -^}f + K 2 E[v\v^\ . (1.12) 
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rV 

Proof: For any x > y, v(x) = v(y) + / v x (z)dz < v(y) + \/x — y\\v x \\2- Now for any a > 0, 

J X 

integrating in y over [x — a,x], av(x) < / v(y)dy H — a 3 ^ 2 1 1 -^a; 1 1 2 ■ Optimizing in a, one obtans 

Jx-a 3 



PI) . Next, 



/ (l + a^Xu-t^fdx < ||u||oo / (l + x 2 )vdx 

J{\x\>C} JUx\>C\ 



< 



< 



VWoa 

C 2 + l 
C 2 1 

C 2 + l 
C 2 1 



'{\x\>C} 



x 2 vdx 



{\x\>C} 



(1.13) 



where we have used (|1.9f) in the last line. Next, since v and v^°°' have the same mass, 

1 



1 

2C 



{\x\<C} 



(v-v^)dx 



1 

2C 



vdx < 



{\x\>C} 



2C 3 



\x\ vdx . 



{\x\>C} 



Hence there exists an xq G [— C, C] such that \v(xq) — t/°°)(2;o)| < 



1 



2C 3 



|x| 2 fdx. Thus, for 



any x G [— C, C], 



{|x|>C} 



3 



< 



2C 3 



E[v\v (co) ] + 2JCE[v\v( c °)] 



where in the last line we have used (11.91) twice. It follows that 



{\x\<C} 



(1 + x 2 )(v - v^fdx < (1 + C 2 ) 



jL(E[v\v^]) 2 + 8C 2 E[vl 



(00)1 



(1.14) 



Combining (jl.lip . (|1.13[) and (|1.14l ). and recalling that M determines C, we obtain the result, and 
see that explicit values of K\ and K2 may be written down. □ 

Hence any result showing that E[v \v (°°)] converges to zero along solutions of (|1.4j) shows that 
v converges to the I • HI 1,1 norm. Such results have been proved in [3] and more recently in 

[7], with an optimal exponential rate. However, in these papers, a weaker form of (|1.9p was used, 
without the second term on the right. Consequently, these papers only deduced the convergence of 
v to t>(°°) in the ||| • ||| 0,1 norm. 

At first sight, it may seem a trivial matter to go from convergence in the ||| • |||o,i norm to 
convergence in the ||| • ||| mj i norm for higher values of m: One might guess that since the steady 
state supported in [— C, C], it should be easy to control the evolution of higher moments 

M2m{v) := L \x\ 2m v{x, t) dx of solutions v of (jl.4p . Upon careful consideration, this turns out 
not to be the case. We shall present a non-trivial argument to show that M^{v{-,t)) is bounded 
uniformly in t in terms of its initial value, but we do not know if this is even true for moments 
higher than the fourth. So while one might expect to be able to prove strong localization estimates 
for solutions of (|1.4p . localization and moment bounds turn out to be somewhat subtle. Using our 
uniform bound on M^v), we shall deduce an exponential rate of convergence of v to in the 
I • ||| m 1 norm for all 1 < m < 2. 
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Moreover, we shall give a new proof of the fact that E[v 1?/°°)] converges to zero at an exponential 
rate. In this new proof, moment estimates play the crucial role: The rate at which .Ej?;!?/ 00 )] 
converges to zero is controlled by the rate at which L |x| 2 vdx converges to / R |s|V°°)(a;)dx. This 
is somewhat remarkable: We are using a simple functional involving no derivatives of v to control 
one that does involve derivatives of v. There turns out to be close interplay be between smoothness 
and localization for solutions of (II. 4ft . and one point of this paper is to explain this interplay, and 
show how it may be used. 

1.2 Convergence of second moments and convergence of the energy 

Let v(x,t) be a solution of (|1.4j) . and let be the stationary solution of the same total mass. 
Then the relative second moment 

a[v\v {oo) ] = a[v(-,t)] -a[v {oo) ] where a[v] = - f x 2 v{x)dx , (1.15) 

2 Jr 

and the relative surface energy 

P[v\v {qo) ] = P[v(-,t)] -f3[v {Qo) ] where p[v] = - f v 2 x dxdx . (1.16) 

2 Jr 

Evidently, the three quantities and /3[v |t>(°°)] are related by 

E[v\vW] = a[v\v^\ + P[v\v^\ . (1.17) 

It follows that if one can show that any two of these converge to zero, so does the third. As we 
have seen above, if one can show that lini£_ ) . 0O E^li/ 00 )] = 0, one concludes as well that 
lim^oo \\v — v(°°}\li y i = 0, from which it certainly follws that lim^oo af^l?/ 00 )] = 0. 

What is perhaps more surprising is that if one can show that lim^oo = 0, one can 

also deduce as a direct consequence that lrm£_ s , 0O £/[v|t>(°°)] = 0, and moreover, one can estimate 
the rate of convergence in the latter limit in terms of the former limit. Let us explain how this 
works, first at the level of formal calculation. 

We easily compute that [1] 

|a[t;(-,t)|t;( 0O )] = -2aK,t)|^ 0O )] + 3/3K,t)|^ 00 )] . (1.18) 
From this and (11.170 we get 

4a[u(-,t)|u (oo) ] = -5a[v{-,t)\v^} + 3E[v(;t)\v^°K (1.19) 
at 

Now, for T > 1 let us integrate both sides of (|1.19p from T — 1 to T to obtain 

a(v(;T)\v^)-a(v(;T-l)\v^) + 5 C a[v(-,t)\v^] dt > 3E[v(-, T)\v {oo) ], (1.20) 

JT-l 

since E[v(-,t)\v^°°^] is monotone decreasing. Now suppose we have an estimate of the form 

\a(v(-,t)\v^)\ <Ke~ xt . (1.21) 
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Using this in (|1.20p yields 

E[v(-,T)\v^] < j(l + e x + je x ^j e~ XT . (1.22) 

In the next subsection, we explain how we shall obtain a rate of convergence estimate for 
a[v(; t)|v (oo) ]. 

1.3 The second Lyapunov functional 

Carrillo and Toscani have made the remarkable discovery [5] that the equation (II, possesses a 
second Lyapunov functional: Define the entropy H[v] by 



H[v] = J R [Y v{x) + 2 \fl v3/2{x) ) 



dx. 



and then the relative entropy by .ff [■uli/ 00 -*] = H[v] — H k/ 00 -*]. We remark that by (jl.llj) . 

E[v] < oo => H[v] < oo . (1.23) 

It is easy to see that v^- 00 ' minimizes H [v] among all non-negative integrable functions v with the 
same mass as v^°°\ and hence H[v \v^°°^] is non-negative. In fact, as shown by Otto |llj . 

H[v\v^]>(J\v-v^\dx\ . (1.24) 

The entropy functional H arises in the theory of the porous medium equation. There is a 
particular slow-diffusion case of the porous medium equation for which the Smyth-Hill densities 
v(°°} are also steady state solutions. This equation, which can be written in the gradient flow form 

is second-order parabolic. For it, the maximum principle applies and provides both positivity and 
uniqueness. Hence H [v [v^ 00 ^] is a Lyapunov functional for the equation fll.25f) . and on account of 
(|1.24p . it is a good one. 

The remarkable discovery of Carrillo and Toscani is that H[v \v (°°'] is also a Lyapunov functional 
for the thin film equation ()1.4|) . and indeed they even show that for strong solutions of (jl.4p and 

T> S, 

rp 

H[v(- ) T)\v {co) } + 2 ( H[v{; t)\v {co) ]dt < F[v(-,S)|v (oo) ] . (1.26) 
Js 

This has the immediate consequence that 

H[v(-,t)\v^] < e- 2t H[v{-,0)\v^] . (1.27) 

The fact that (jl.26p holds for solutions v of the slow diffusion equation (jl.25p was discovered 
by Otto [TT] using methods from the theory of optimal mass transportation, and in particular, the 
notion of displacement convexity, as we shall explain in more detail in the next section. 
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The fact that (jl.26p also holds for solutions of (|1.4p is far from obvious, but was proved by 
Carrillo and Toscani using the fact that the equation (jl.4p can be written as 



vt 



y + Kv) 



+ v 



y + Kv) 



with h(y) = ^/Qv 1 / 2 and <J>(i>) = vh'(v). Combining (|1.27p and (|1.24p . they then deduced 



\v - u (oo) |dx < e~ < -y/if[t;(-,0)|u( 00 )] , 



(1.28) 



and raised the question of proving that the energy E[v(-,t)\v^°°^] converges to zero. 

As explained in the previous subsection, to do this, it suffices to prove that 
converges to zero. We shall do this by proving in the second section of this paper an inequality 
relating for the entropy iJfulin 00 )] and a[u(-, £)|u(°°)]. We shall prove: 

1.4 LEMMA. For any non-negative integrable function v onM with a finite second moment, let 
■y(°°) be the Smyth-Hill density with the same mass. Then 



a[v\v {oo) ] < 2^ a[v(°°)]y H[v\v(°°)] + H[ 
Granted this lemma, the bound (|1.27p now gives us 



\v\v 



(oo)i 



a 



V(;t)\v^] < e"*2y 'a[v(°°)]y/ 'h[v(;0)\v(°°)] 



+ e- 2t H\v( 



>(°°)l 



and then from (jl.22p we have E[v(-,t)\v^°°^] < Ce~* for an explicit constant C. Finally, by 
Lemma II. 3\ we obtain exponential convergence in the ||| • norm. This outlines the general 
scheme of our strategy. 

We remark that in [7] it is shown that for a certain class of solutions, E[v(-,t)\v^°°^] decays like 
e~ 2t , which is twice the rate implied by our result that £)|v(°°)] decays like e~~ t . However, our 
construction is somewhat less delicate, and the information we obtain on moments and localization 
is new. 



1.4 Properly dissipative weak solutions 

The theory of the thin film equation is not yet in a well developed state. Basic issues of existence 
of strong solutions and uniqueness remain open. For which class of solutions can the formal cal- 
culations above be made precise and rigorous? We now introduce such a class of weak solutions, 
called properly dissipative weak solutions. 

The key estimates used in our convergence analysis depend on the fact that a[i>(-, t)\v (°°'] 
satisfies (fl~19]) and that H[v(-, t)|v (oo) ] satisfies ([06]) . Therefore, what we need is existence of 
weak solutions with these properties. 

1.5 DEFINITION. A non-negative measurable function v 6 C°(M x [0,oo)) such that for some 
fixed M, the mass, and each t > 0, 



j v(x,t)dx = M and E[v(-,t)] < oo , 
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is called a weak solution if for all £ £ Cg°(M x (0, oo)) 

/ / (-v( t - 3(v x ) 2 ( xx - 2vv x ( xxx ) dtdx = 0. (1.29) 

J JRx(0,oo) 

It is called a properly dissipative weak solution if moreover: 

• For all * > 

2/ {{(vlv^dt + Hlvi^t)^} < H[v(-,0)\v^}, (1.30) 
Jo 

• For all t > 1, 

a[v(-, t)\v {Qo) ] - a[v(-,t - l)\v {co) ] + 5 f* a[v(-, s)\v {co) ]ds > 3E[v(t)\v {oo) ] . (1.31) 

Our main existence theorem for properly dissipative weak solutions is the following: 

1.6 THEOREM. Let v$ be any non-negative integrable function on M such that E[vq] < 00 and 
such that 

-^4(^0) = / x 4 vo(x)dx < 00 . 



Let M be the total mass o/vq, and let denote the Smyth-Hill steady state with the same mass 
M. Then there exists a properly dissipative weak solution v such that v(x,0) = vo(x). Moreover, 
there is an explicit constant K3 depending only on M and E[vq] 

[ \x\ 4 v(x,t)dx < [ \x\ 4 v (x)dx + KzHlvolv^] . (1.32) 

As noted in in (|1.23p . i?[uo|u^ 00 ''] is finite whenever E[vq] is finite. Theorem 11.61 is proved in the 
third section of the paper. We now state another of our main results: 

1.7 THEOREM. Let vq be any non-negative integrable function on M such that E(vq) < 00 
and such that M^vq) < 00. Let M be the total mass of vq, and let i/ 00 ) denote the Smyth-Hill 
steady state with the same mass M . Then for any properly dissipative weak solution of \l-4\ ) with 
v(x,0) = v (x), 

lib - ^ (oo) lll?,i < K A e- 1 (1.33) 
where K4 is a positve constant depending only on M and E(vq). Moreover, for all 1 < p < 2, 

\\v- v (oo) ||p,i < ^ 5 e- (2 - p)t (1.34) 
where K§ is a positive constants depending only on p, M, M^vq) and E(v$). 
This theorem is also proved in the third section. 
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2 Mass transportation and the thin film equation 

We rely in an essential way on methods of optimal mass transportation to both construct and 
analyze our weak solutions. In this section we briefly recall the points that are essential here. See 
|13j or m for more information. 

For M > 0, let Ai m denote the set of non-negative Borel measure [i on R with fi(M) = M and 
such that 

/ |x| 2 d^(x) < oo . 
Jr 

For fi, v G M-m, define v) to be the set of all non-negative Borel measures 7 on R 2 such that 
for all Borel sets icM, 

7(A x R) = n(A) and y(R x A) = v(A) . 

The set T(fi, v) is the set of all couplings of /i and v. 

The 2-Wasserstein distance between [jl, u G Mm, W2(/i,f), is defined by 

W 2 (^) = inf (/ |x-y| 2 d 7 (x,y)) . (2.1) 

By the Brenier-McCann Theorem, when fi and v are absolutely continuous, the infimum is 
attained at a unique optimal coupling 7*, which is concentrated on the graph of the derivative <p x 
of a convex function <p: For all Borel sets 4,5cR, 

7 (ixB) = ^n^(B)) . (2.2) 

It follows that 

W 2 Gu,^)= / \x - <p x (x)\ 2 dn(x). (2.3) 

We write ~W\{v,w) to denote W^vdx, wdx), as is standard. 

The fact that the optimal coupling is induced by a map x 1— > <p x (x), which is ^-almost everywhere 
invertible, yields an interpolation between fi and v. Let ip be the a convex function such that (|2.2p 
defines an optimal coupling of \i and v. For t G [0, 1], define the convex function y>W by 

2 



<pV( x ) = (i-t)^- + Mx) . 



Then ip$ interpolates between the identity and ip x , and we define and G M 



M 



/' 



it) 



{a)=^)-\a)). 



The map t 1— )■ //W i s McCann's displacement interpolation between and ^, 

There is another way of expressing this that will be useful to us. Let T : R — > R be measurable. 
Then Tftfj,, the push-forward of \x under T, is the measure given by (T#fi)(A) = fi(T~ 1 (A)). Thus, 

^ = (V W )x#^- 

A functional Q on Mm is sa -id to be A displacement convex, [3 [13] if for all fx, v G A4jv/, the 
displacement interpolation of /i and f satisfies 

A*(i - t)w 2 (/x, 1/) + gQjfi)) < (1 - t)g(|i) + tg(u), (2.4) 
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for all < t < 1. In this case, 

AW|(m, i/) + Inn <Q[ V ]- gfa] . ( 2 .5) 

Then, if fj, minimizes Q so that the subgradient of Q vanishes at \i (seepQ), (|2.5p reduces to 

AW| (//,*/) < - g\p\ , (2.6) 

which is known as a Talagrand inequality for the functional £/. 

This is of immediate relevance to the entropy functional H[v], regarded as a functional on Mm 
in the obvious way, since this is 1-displacement convex as discussed in [5], and is minimized by 
v(°°^dx. Thus taking d/i = v^°°^(x)dx and du = v(x)dx in Mm-, (|2.6p specializes to 

V%M<H[v\vW] . (2.7) 

Proof of Lemma 11.4b Because of (12, 7p . it remains to show that 

a[v\v {oo) ] < 2^/a[ V (°°)]W 2 (i;,t; (oo) ) + W|(w,v (oo) ) . 

To see this, let 7 be an optimal coupling of v(x)dx and t> (°°\x)dx, and note that ^/ a[v] = 
1 1 ^ 1 1 z, 2 (M 2 ,cd-y) an d \/ a[v(°°)] = ||y||L 2 (R 2 ,d7) where we write x and y to denote the functions (x, y) 1— > x 
and (x, y) 1— >• y respectively. Then by the triangle inequality, 

\y/^v\-y/a[ v (<x>)]\ < \\x-y\\ L 2 [R 2 :dy) =W 2 (v,v^) . (2.8) 

□ 



3 Construction of properly dissipative weak solutions 

Define /Cm to be the subset of Mn consisting of absolutely continuous measures v(x)dx such that 
E[v] < 00. 

We now implement the JKO scheme [TJ [6]: Given a small time step r > 0, and given an initial 
datum vo £ /Cm, we recursively define a sequence (u^)neN by Vq = vq and 

< +1 = argmin,^ jr^M^ 00 )] + , (3.1) 

The first step of our analysis is to show that the variational scheme in (|3.ip has a solution and 
to derive the Euler-Lagrange equation for this variational scheme. 

3.1 LEMMA. Given r > 0, let vq £ /C . 

(i) The discrete variational scheme \3.1\) admits a solution {v£}neN- 

(ii) Each satisfies 

K) xxx £ L£ C (R), (3.2) 

and 

vl(a) = «)x(a) = 0, a £ R. (3.3) 
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(Hi) Let ijj' n be the the optimal transportation plan such that vj l _ 1 (x)dx = (ijj' n )^j=vj l (x)dx. Then 

tp' n (x) =x + t [xv T n {x) - v T n {x){v T n ) xxx {x)\ . (3.4) 

The proof of this lemma is closely patterened on a proof of Otto |10| . We present some details 
for the convenience of the reader since Otto's energy functional differs from our in having the term 
(1/2) J R x 2 v(x)dx replaced by the Lebesgue measure of the set {x v(x) > 0}. The reader familiar 
with |10j . or even [6], may wish to skip ahead to Lemma 13.31 For other readers, we point out that 
(|3.4p is the Euler-Lagrange equation for the discrete time variational principle, and through (|3,4p . 
one makes a direct connection with the thin film equation (jl.4p . 

Proof of Lemma I3.lt Proof of (i): It is enough to show that for given vq G ICm, there exists a 
solution of 

v\ = argmin^^ lrE[v\v^\ + \w$(vo,v)\ . (3.5) 

Let v k ke~N be a minimizing sequence in (|3.5p . As {E(v k \v ^°°^)}kef>i 1S bounded, there exists 
v% : R — > [0, oo) such that 

v k — > vi, locally uniformly on K (3.6) 

for a subsequence and that 

£(^i|v (oo) ) < Hminf E(v k \v^). (3.7) 

v(xi) - v(x 2 ) 



Indeed, by (jl.lip and the elementary bound sup — — < y2E\v\, the Arzela- 

xi,x 2 eM.,xi^X2 V \xi — x 2 \ 

Ascoli theorem may be used to prove (|3.6j) . Due to (|3.6|) and the boundedness of {(i ,fc )'}feeN m 

L 2 (R) we have v[ G L 2 (R), (v k ) x (vi) x , in L 2 (R), and in particular 



[ \(vi) x \ 2 dx< hminf / \(v k ) x \ 2 dx. (3i 



Using the inequality (|2.8p . for all vo,v\ G /C, we deduce that the boundedness of {W 2 {vq, v fc )}fceN 
implies that {f R x 2 v k dx} k N is bounded. Then, by this, (|3.6p and Fatou's lemma, we deduce 



/ 



x 2 v\{x) dx < oo. 



Also, f^ r v k (x)dx = M - f R \i_ rjr] v k dx > M - ^ J R x 2 v k dx and so lim rtoo J^ r v k dx = M, 
uniformly in k. And thus, J R v\dx = 1. Combining, we have E(vi\v^°°^) < liminffc_ ! , 00 E(v k \v^°°^) 
along a subsequence. Next, the fact that 

Wi(v , Vl ) < liminf Wi(v , v k ). (3.9) 
k~ >oo 



follows from the lower semicontinuity of the Wasserstein distance, see [10] or [13j . This concludes 
the proof of (i). 

Proof of (ii) and (Hi) The key step here is a variational argument of Otto [10] showing that, 
in our case, at each n the optimal transportation plan ijj' n such that v 1 n _ 1 {x)dx = (ip' n )^v^ l (x)dx is 
related to v T n through 

1 

T 



(x - ip' n (x))£(x)vl(x) dx + 



3 

x^{x)v r n {x) - -«)^'(i) - v T n {v T n ) x i"{x) 



dx = 0. (3.10) 
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To make the variation of v^, consider any £ £ Cq°(R), and define 7j(a;)da; = (Id+e^)^(v^ l (x)dx). 
We now work out the effects of this variation on each term in the functional being minimized in 
(13, 5p . starting with the Wasserstein distance. 

Let 7 denote the optimal coupling of vl n {x)dx and vj l _ 1 (x)dx. Then 7 := ((Id + e£) ® 7<i)#7 is 
some (non-optimal) coupling of vdx and u^_ 1 (x)dx, and hence 



1, 



1 



W 2 2 (6,<-i) " ^ e / ^-^(x))£(xKdx + 0(e 2 ) 



2 z v . 2 
As for the second moments, we have 



- / x 2 v{x)dx — — J x 2 v 1 n {x)dx 



x£(x)^(a;)dx-)-0(e 2 ). 



Next, since 



d.r 



\ f vldx-- f {v T n ) 2 x dx 

1 JR 1 JR 



dy 

■£ 



dx 



with ^. = i-eC'(y) + 0(e 2 ) 



dx + 0{e A 



(3.11) 

(3.12) 

(3.13) 
(3.14) 



Now (l3TTUj) follows by combining (pTTT]l . (l3TT2l . (^T4]) and replacing £ by 

Next, we complete the proof of (ii): As in Proposition 1.3 of \W\, because any v £ /C is 
continuous, the set {u^ > 0} consists of countably many intervals. Let / = (a, 6) be such an 
interval. Then we have from (|3.1U|) that for some / £ L^ c (R) (expressible n terms of ip' n ), whenever 
£ that is smooth and compactly supported in /, 

"3 



l{vlfJ{x)+vl{vlU"{x) 



dx 



£{x)f(x)v T n {x)dx 



Hence, taking distributional derivatives, 



{Ox - 



v T J, in D*(I) 



(3.15) 



It follows that v T n £ Hf^(I), see [10] for details. For any v £ H^ c (I), 

*$VxVxx and i^^x^xx ~~ ^^x^xx ~t~ Wxxx 



" ( \2 
nK v x)x 



(3.16) 
(3.17) 



Hence we can rewrite (|3.15p as 

(Vn)xxx(x) = f(x) On I , 

for the same /. Next, as noted above, {v^ > 0} is a union of countably many intervals {(cij, bj)}j>i 
Integrating by parts, using the fact that v T n vanishes on the endpoints of each interval, 

3 



&nf dx 



K)ie-vi(viu" 



dx 



E 

00 

E 



dx 



Vn(Vn)xxx£,dx 



1 h C°i 

.7=1 ^ 



(3.18) 
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Equivalently, 

1 oo 



2 

i=i 



Since £ is an arbitrary smooth compactly supported function, this implies that (v^) x (aj) = 
(vl) x (bj), for all j > 1. 

Finally, we have enough regularity to deduce (Hi) from (|3,10p . which implies that 



ijj' n (x) = X + T 



3 

xv T n {x) + -«)L0) - ««)*)a«(a?) 



Simplifying this using (|3.16p . we obtain (|3.4p . □ 
At this stage, we depart from Otto's analysis in [10]. Our next goal is to show that the entropy 
iTfi^l?/ 00 -'] is monotone decreasing in n and to relate the decrease to what what one would guess by 
a formal differentiation argument with the continuous time evolution equation. For this we exploit 
the displacement convexity of the functional H[v \v (°°)]. First, we make a definition: 

3.2 DEFINITION (Entropy dissipation). The functional D is defined on TC M by 

D[v] := / (x+ V6v x (x) I v(x) dx . (3.19) 

J Ml 

3.3 LEMMA. Fix n G N and let be a solution of §3. Then, the following inequality holds 

< -D[v n ] - — \ (v n ) {v n ) x 

T 24 J R 

< -2H[v^\v {oo) ] (3.20) 

Proof: The second inequality follows from the first using the entropy-entropy dissipation inequality 
2H[v\v i - cc ' S} ] < D[v], as explained in [5]. Hence we must prove the first inequality. 

Let < t < 1, denote the displacement interpolation between v T n and v T n _ v Since H[v] 

is a displacement convex functional, it follows that for t € (0, 1) 

H[vl_ x ] - H[v T n ] > \ (H[(vZ)®] - HK\) , 

and moreover, the right hand side is monotone decreasing as t tends to zero . 
By a standard computation, [T| I13j. 



Jim \ " HKW^i) > ~ I (f [<]) («. 



x) T dx 



where 

x 



5H ~ 2 



Integrating by parts and using the Euler-Lagrange equation (13. 4h . 

K_^°°] - #K> (oo) ] > /[^ + v / 6 v / <],«(rE)-xKdx 

./TCP ^ 



/ [x + v / 6( A /<)a : ](2; - {v T n ) xxx )v T n dx, (3.21) 
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where we have used (|3.4p . Adding and subtracting x + VQ(y/v^) x , we deduce that 



HK-i\v (co) ] ~ HK\v {oo) ] >r \x + V6(J^) X 



where 



J := -r 



vL dx + J 



vZ dx. 



(3.22) 



(3.23) 



Integrating by parts, using Lemma [3.41 proved just below, to justify certain of these integration by 
parts, we obtain 



J=^r[(v^dx + [K)^K)l x dx 

"J JR z JR 



+ r 

24 



K)-v 2 K) x dx 



each term of which is non-negative. Combining (|3.24p and (|3.22p we obtain the result. 
3.4 LEMMA. Let be the nth step in a solution of \3. 1}) . Then, 



K) 1/s 



g C(R), 



2/ " n(a) = 



(3.24) 
□ 

(3.25) 



Proof: Without loss of generality assume that v^(0) = = (t>^) x (0). Assume on the contrary that 
> <5 3 > for some 5 > 0. This implies (t^(x)) 1 / 2 > Cix 3 / 5 for some constant C\ > 0. Indeed, 

we easily deduce from the assumption that ^j§fj > <5 > which implies easily that [(f^) 5 ^ 6 ]a; > C5. 
Then, (t;^(x)) 5 / 6 > C5x with C > holds. This immediately implies the assertion. Now, by (13.2p 



we get («£)a; < C 2 x, for some constant C 2 > 0. Hence, < C^/s = Cx 12 / 5 -)• as x ->■ 0. 



This is a contradiction. 

We next examine the behavior of second moments along the discrete scheme. 



□ 



3.5 LEMMA. Fix N, M G N with N > M and let v T n be a solution of [3l\) . Then, the following 
inequality holds true. 



N 

a[v T N \v {oo) } - a[v T M \v (oo) ] + 5r a[v]\v^} 

j=M+l 
N 

> 3r E[v]\v^}-tE[v \v {co) } . 

j=M+l 

Proof: Since 

aK\v {oo) ] - aK_!k (cx>) ] = aK] - a[v T n ^] = [ (\x\ 2 - |<(x)| 2 K(x)dx , 



(3.26) 



a simple computation using the Euler-Lagrange equation (I3.4p and the regularity results to integrate 
by parts, one obtains 



a\vZ\vW] 



a v r 



-il^ (oo) ] > -2r«K_ 1 |t;( 00 )] + 3r/3K_ 1 |t- 00 ] - JwlKX-l)- 



(3.27) 
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Thus, 

TV TV 1 TV 



a[v- N \v^]-a[vl I \v^] + 5r £ a[vj\v^]>3r E[v]\v^]-\ £ Wjiv^v]), 

j=M+l j=M+l j=M+l 

(3.28) 

which is obtained by summing up the estimate (|3.27D . 

Now, observe that thanks to the variational structure of (|3.ip . we obtain the following estimate 
for free 

TV 

E W 2 2 (^-i>^ T ) ^ 2t(e[v \v^} - E[vjf\v^]) t (3.29) 

k=l 

and hence using this in (13.28H we obtain (|3.26p . □ 
We next control the fourth moments. 

3.6 LEMMA. Fix n € N and let be a solution of $3.l\) with initial data vq € /C having a finite 
fourth moment; i.e., M^{vq) < oo. Then M 4 (t^) is bounded uniformly in n. Indeed, there is an 
explicit constant K% depending only on M and E[vq] such that 

M 4 «) < M A (y Q ) + K 3 H[v \v^] . (3.30) 

Proof: By the displacement convexity of the functional M 4 (/) we deduce that 

M 4 «_ 1 ) > M 4 «) + ^- [(x + r(x-(vl) xxx )) 4 vldx 

dT JR r=0 



= M 4 «) + 4r / x 3 (x - {v T n ) xxx ))v T n dx 

JR 

= M 4 «) + 4rM 4 «) -t x 3 {v T n ) xxx v T n dx . 

JR 

(3.31) 

We now integrate by part on the last term, using Lemma [3. 41 to justify the calculations. We obtain 
M4(<) ~ M4(< - l} <-4M 4 «)-l2 / K) 2 dx + l8 [ x 2 [(vl) x ] 2 dx. (3.32) 

T JR JR 

The last term on the right hand side of ()3.32|) is estimated as follows: 

x\vi) x fdx = [ x'KfKr^iKu'dx 

JR 

< ii<ii^ 4 (M4«)) i/2 (^«)- 3/2 «)^xy 2 , 

< 4M 4 «) + M|*L / (v^ 2 (vl)idx . 

1,3 JR 

Combining this with (|3.32p and using (jl.lip and Lemma 13.31 we obtain 

M 4 «) - M 4 «_0 < K,{HK^\v^\ - HK\v^]) . 
Telescoping the sums gives the result. □ 
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We are now ready to prove Theorem 11.61 the existence of properly dissipative weak solutions: 

Proof of Theorem II. 6t Define v^ T \x,t) = v^(x) for nr < t < (n + l)r. It is then standard to 
show [H I10| [6] that from the family {f^ T ^} T >o, one can extract a weakly convergent subsequence, 
and that the weak limit is a weak solution of (jl.4p in the sense of f|1.29|) . See [10] for such an 
argument. 

Next, for any r > and t = nr, we have from Lemmg [3T3l that 

2 f l H(v\v {co) )dt + H[v{-,t)\v {oo) ] < H[v{-,0)\v (oo) ] . 
Jo 

A standard convexity and lower semicontinuity argument shows that this inequality is preserved 
along weakly converging subsequences. This proves (|1.3Up . The argument is relatively straightfor- 
ward since the "small" side of the inequality is a weakly lower semicontinuous function of v(x, t), 
while the "large" side of the inequality only depends on the initial data. 

A more involved argument is required to prove (jl.3ip since in this case the solution v(x,t), and 
not only the initial data, occurs on both sides of the inequality. We therefore need continuity, and 
not only lower semicontinuity, of the functionals on the "large" side. 

This is provided by the uniformly bounded fourth moment. On account of this, we conclude 
that along any weakly convergent subsequence, {v^^} 

lim a[v {Tk \-,t)] = a[v(-,t)] , 

k— too 

while again a standard convexity and lower semicontinuity argument shows that 

lim E[v {Tk \-,t)] > E[v(-,t)] , 

k— ¥oo 

Since E is decreasing along each solution of the discrete scheme, (jl.3ip holds with i>( Tfc ) in place 
of v , and then by what we have said above, the inequality is preserved in the limit. This proves 
(TL3T|) . □. 
Proof of Theorem 11.7): Since v is a properly dissipative weak solution is satisfies (|l,30p . and 
hence 

H[v(-,t)\v^} < e- 2t H[w\v^} . 

Then by Lemma ll.4| 



a v 



where K depends on v$ only through M and H[vq\v (°°^]. Then since v is a properly dissipative 
weak solution is satisfies (|1.3ip . and hence v satsfies (|1.22p . This proves (|1.33|) . 
Next, since \x\ 2p = (|x| 2 ) 2_p (|x| 4 ) p_1 , Holder's inequality yields 

\x\ 2p \v-v {oo) \ 2 dx< I I \x\ 2 \v- v {co) \ 2 dx) P ( ( \x\ 4 \v-v {co) \ 2 dx XP 



By the pointwise bound, there is a constant C depending only on E[vq] such that \\v — ||oo < C, 
and hence for another constant C depending only on E[vq] and M^Do], L |x| 4 |t> — v^°°'\ 2 dx < C . 
Now (fL34"l) follows easily. □ 
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